Cooling a micro-mechanical resonator to its ground state by measurement back-action 
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We present an analysis of the cooling of a micro-mechanical resonator by means of measurements 
and back action. The measurements are performed via the coupling to a Cooper-pair box, and 
although the coupling does not lead to net cooling, the extraction of information and hence entropy 
from the system leads to a pure quantum state. Under suitable circumstances, the states become 
very close to coherent states, conditioned on the measurement record, and can hence be displaced 
to the oscillator ground state. 



I. INTRODUCTION 

Micro-mechanical resonators are elastic micron-sized 
beams fixed in one or both ends with tranversal modes 
of vibration with frequencies in the range 10 to 1000 
Mhz [l[. Applications of these oscillators include ultra- 
sensitive force detectors, electrometers and displacement 
detectors. Recently, interest has been directed towards 
developing schemes to make these resonators display non- 
classical behaviour. Here, one tries to cool the resonators 
to their vibrational ground-state 0, [H, H land to generate 
different non-classical states of motion l|, || . The quan- 
tum regime represents the fundamental resolution limit 
for these systems when they are used for ultra-sensitive 
detection purposes, and by making special, squeezed 
states or macroscopic superposition states, one may sig- 
nificantly enhance their performance. When operated in 
the quantum regime, mechanical resonators may find ap- 
plications as quantum information storage devices and 
as mediators of quantum states between other quantum 
systems. Also, proposals exist to produce quantum corre- 
lated states of motion of two oscillators [6[ and to estab- 
lish quantum entangled states of hybrid systems involv- 
ing a mechanical resonator and, e.g., a quantized field Q 
or an atomic cloud Q. 

For all of these ideas to be realized, however, it is neces- 
sary to coherently control the quantum state of these sys- 
tems, and with their low vibrational frequencies, even at 
cryogenic temperatures, their thermal excitations are sig- 
nificant, and the quantum state will in general be a ther- 
mal mixture. Different means of cooling have thus been 
proposed and in particular mechanical resonators which 
also act as optical resonators, and which can be detuned 
with respect to the frequency of incident laser radiation, 
have been cooled in analogy with frictional and sideband 
laser cooling of atoms and ions 0, [l(| [H| ■ Other means 
for cooling, e.g., by the back action of the coupling to 
resonant circuit elements [l2[ have been analyzed and 
demonstrated, and in this article we present a way to 
turn a mixed thermal state into a pure coherent state of 
motion using the quantum mechanical back-action due to 
measurements on the mechanical oscillator. We imagine 
these measurements to take place by a Cooper-pair box 



which is capacitively coupled to the resonator and which 
is read out by a single electron transistor (SET). The 
measurement does not cool the system, but it extracts en- 
tropy and provides a pure coherent state, and hereafter, 
or along with the meaurement process, a simple feedback 
mechanism may be used to displace the state, so that the 
mechanical ground state is eventually obtained. In our 
simulations we are able to estimate the efficiency of this 
method and to incorporate various effects such as finite 
interaction time and coupling to a thermal reservoir, and 
see how they affect the cooling procedure. In Sec. II, we 
introduce our physical system. In Sec. Ill, we describe 
the measurement protocol and the formalism describing 
the effect of measurements on the oscillator. In Sec. IV, 
we present our simulation results. In Sec. V we describe 
the use of feedback to reach the oscillator ground state, 
and we show how well the measurements and feedback 
are able to balance the heating of the oscillator due to its 
coupling to a thermal reservoir. Sec. VI concludes the 
paper. 



II. THE SYSTEM 

The micro-mechanical resonator studied is a cantilever 
beam. This is a beam rigidly attached in one end and 
completely free in the other end. The beam has flexural 
modes of vibration in the directions perpendicular to the 
axis of the beam. The lowest frequency mode can be 
modelled as a harmonic oscillator and one can assume 
weak coupling between this mode and the other modes 
of the cantilever [l[. 

The other part of the system is a Cooper-pair box. 
This is, effectively, a charged capacitor with an anhar- 
monic ladder of discrete charge states, and where a two- 
level system can be effectively identified with two specific 
charge states, |— ) = \n) and |+) = \n + 1). These states 
correspond to n or n+1 Cooper-pairs being present on the 
capacitor "island" of the Cooper-pair box, respectively. 
When restricted to these two states, the Hamiltonian of 
the Cooper-pair box can be written [l3l |: 



Hcpb = E c 5na z —a x 



(1) 
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FIG. 1: This figure shows a sketch of the setup. The double 
arrow shows the direction in which the oscillations are imag- 
ined to take place. V™ and Vg are the voltages applied to 
the micro-mechanical oscillator and the other gate electrode, 
respectively. C™ and Cg are capacitances between the elec- 
trodes and the Cooper-pair box "island". 



arbitrary initial superposition state = c±\—) + c 2 |+) 
we transfer the system to the new superposition state 
(ci cos(<5) + c 2 i sm(S))\-) + (c 2 cos(5) + ci£sin(<5))|+). By 
switching on an off Ej for short intervals of time, and 
by letting the Cooper-pair box and cantilever oscillator 
interact, the two systems become entangled. Second, us- 
ing a radio-frequency single electron transistor (rf-SET), 
it is possible with high sensitivity to measure the charge 
state of the Cooper-pair box [l[ . Currently, with realistic 
Cooper-pair box parameters it is possible to measure the 
charge of the Cooper-pair box within a time r m = 4ns 
[l6|. This is much shorter than the cantilever time scale 
of motion, and we will for simplicity of our analysis as- 
sume that measurements are instantaneous with respect 
to the dynamics of the oscillator. Recent experiments 
have demonstrated an alternative read out of a Cooper- 
pair box by its modification of the transmission proper- 
ties of a stripline cavity (l7j |. By carrying out measure- 
ments of the charge states of the Cooper pair box, we can 
extract information, and thus remove uncertainty about 
the cantilever motion. 



Here Eq is the so-called charging energy of the Cooper- 
pair box, Sn is a tunable dimcnsionless parameter taking 
values between -1/2 and 1/2 and Ej is the Josephson 
energy. The operators a x , a y and a z are defined by their 
action on the states |— ) = |n) and |+) = \n + 1) the 
same way as the corresponding Pauli operators acts on 
the eigenstates of the projection of a spin onto the z-axis. 
If Cooper-pair box parameters are chosen [lj such that 
the second term in the Hamiltonian can be neglected, the 
charge states are eigenstates of the system. 

In Fig. Q]a schematic sketch of the setup is shown. We 
see that by putting an electric voltage on the cantilever 
we will get a capacitive coupling between the systems. 
The coupling between the Cooper-pair box and all other 
modes but one of the two lowest frequency modes will 
be weak and can thus be neglected [lj]- The coupling 
between all modes is also weak. Our system therefore 
consists of a Cooper-pair box coupled to a harmonic os- 
cillator. When Ej = 0, the Hamiltonian of the combined 
system is [1(: 

H = fiwa) a + E c 5na z + A(a^ + a)a z (2) 

where w is the frequency of the oscillator and where A 
is the coupling strength between the systems, depending 
on the displacements of the cantilever by x oc a + a^ from 
its equillibrium position and being proportional to the 
voltage between the systems [l|. 

The Cooper-pair box has two advantageous properties. 
First, Ej can be switched to values much larger than 
E c Sn for a short period of time At [l^]. The box is 
then said to be taken to its degeneracy point and coher- 
ent transitions are driven between the charge states. For 
E 2fr l equal to the angle 5, we say that we have applied 
a J-pulse causing a Rabi oscillation between the charge 
states of the Cooper-pair box. This is a unitary process, 
and applying the 5-pulse to the Cooper-pair box in the 



III. MEASUREMENTS 

We want to gain information about the state of the os- 
cillator, and we propose the following measurement pro- 
cedure. At first the Cooper-pair box is in either one of 
the charge states. Then a fast 7r/2-pulse is applied to the 
Cooper-pair box transferring it almost instantaneously 
into an even super-position of the charge states. After a 
time interval r during which the Cooper-pair box and the 
cantilever interact, a second 7r/2-pulse is applied to the 
Cooper-pair box, and then the charge of the Cooper-pair 
box is measured, returning it to one of the charge states, 
so that the procedure can be repeated over an over again. 

We can give the following picture of how the mea- 
surement procedure provides the information about the 
state of the oscillator: After the first 7r/2-pulse puts 
the Cooper-pair box in a super-position of the charge 
states, the interaction Hamiltonian terms containing a z 
will build up a relative phase A<f> between the two charge 
states. Since one of the terms containing a z in the Hamil- 
tonian depends on x oc (a* + a) the phase Acj) will 
depend on the state of the oscillator mode. The second 
7r/2-pulse will rotate the state of the Cooper-pair box so 
that the relative phase A(j) is converted into a population 
difference, and hence, a measurement of the charge state 
populations yields the information about the oscillator 
state. To properly account for this information retriveal 
and its consequences for the dynamics of the cantilever 
quantum state conditioned on the read-out sequence, we 
describe the dynamics and the measurements in a den- 
sity operator language. Let p be the density operator of 
the motional oscillator state prior to a measurement, and 
assume that the Cooper-pair box is in the |— ) state (the 
treatment of the Cooper-pair box being initially in the 
|+) state is fully equivalent). After the first 7r/2-pulse 
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we have for the combined system the density operator: 



P 



P®(|-)(-|+i|+X-|-i|->(+| + |+)<+|) 



(3) 



After an interaction time r the density operator becomes 
UpU^, where U — exp( ~ l ^ T ) with the Hamiltonian in 
(J2j) . Since the charge states are eigenstates of (J2]) we 
can get a simple expression for U pU^ by introducing the 



operators U± = exp(— iojr^a ± ( 



E c Sn 



+ nla 1 + 



&)))), 



acting only on the oscillator space. Using this notation, 
we get U(jK»Qa)(J3\))W = l7 Q pC^®(|o)(/3|), with a,p = 
+ and — , and UpU' can hence be written 
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UpU^ = 

>_®(|-X-|)+v+-®(l+X-l) 

-tp-+®(|-X+l)+P++®(l+X+l)) 



(4) 



with p a f3 = U a pUl. After the final 7r/2-pulse on the 
Cooper-pair box we have the state: 



~i{p-- 
+i(p— 
(P- 



VP- 
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- P+- 

- P+- 
P+- 



- P+H 
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)®(|-)(-|) 

)®(|-)(+l) 
)®(l+)(-|) 

® (l+X+D) 



(5) 



At this stage the projective measurement is performed 
on the Cooper-pair box charge state basis, i.e., the state 
is replaced by the first (last) line of this expression if the 
Cooper-pair box is found in state | — (+)). Introducing 



operators Mz, 



U-TU+ 



the effective modification of the 



cantilever state due to the interaction with the Cooper- 
pair box and its subsequent detection can therefore be 



written p — > p^ 



FT 



occurring with the proba- 



bilities Pzr 



M T p)' 

ti(M^-M T p). The operators U± are expo- 



nential in a, a* , and a^a, i.e., they are analytically known 
combinations of a displacement and a rotation in the os- 
cillator phase space and can be readily evaluated, e.g., in 
the eigenstate basis of the oscillator. This dynamics can 
thus be readily simulated on a computer. 



IV. SIMULATIONS 

Having derived the measurement operators Mzp we can 
numerically simulate the measurement procedure using 
random numbers to decide between the measurement 
outputs in each step of the measurement sequence de- 
scribed above. 

Initially, the system is assumed to be in thermal equi- 
llibrium with a heat reservoir of temperature T. The 
energy difference between the Cooper-pair box charge 
states is assumed to be large compared to kT. Hence, 
we can assume that it is initially in the state — ) and 



the initial density operator can therefore be written 

p = p® (|-)(-|) with : 



E 



exp(hiu/kT) — 1 
exp(hcu I kT) 



e^q>(—hojn/kT)\n)(n\ (6) 



In the following we take an oscillator frequency of u> — 
50MHz and an initial temperature of 50 mK, so that the 
oscillator occupies the lowest hundred oscillator states. 
In our simulations we assume an interaction time r = 
0.1/u) and the charge state energy separation is given by 
E c 5n = IQQhuj. In our simulations we determine the con- 
ditioned dynamics of the cantilever over 160 subsequent 
interaction and measurement sequences, and we repeat 
such simulations 200 times to be able to characterize the 
typical properties of our detection records. Note, how- 
ever, that we do not form an ensemble averaged density 
operator by simply adding together the individual simu- 
lation records: the experiment is supposed to take place 
on a single cantilever, and the specific detection record of 
a single run is assumed to be available to the experimen- 
talist, who can thus infer the cantilever quantum state 
according to the above analysis. 

We first address the approach of our cantilever to- 
wards a pure quantum state as a consequence of more 
and more information about the state becoming known 
via the measurements. For this purpose we calculate in 
each run the von-Neumann entropy of the time depen- 
dent density operator, and in Fig. [2] we plot the average 
over all realizations of this quantity as function of time. 
As is clearly seen, the mean value of the von-Neumann 
entropy decreases with the number of measurements per- 
formed and tends to zero after sufficiently many mea- 
surements. After a hundred measurements the entropy 
is practically zero and we have full knowledge of the pure 
state of the oscillator in every realization. 



& 2.5 - 
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FIG. 2: The averaged von-Neumann entropy of the can- 
tilever state conditioned on the measurement record is plotted 
against the number of measurements performed in over 200 
independent simulations. Parameters used are uj = 50MHz, 
A = 1.5hcu, E c Sn = lOOfiw and an initial temperature of 
T = 50mK. 
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Having thus asserted the successful convergence of the 
quantum state towards a pure state, i.e., a wave func- 
tion, we will now identify what state is actually pro- 
duced. It is natural to quantify the state by the posi- 
tion and momentum variables, and in Fig. [3] we show 
the uncertainty product of these observables, averaged 
over 2000 simulations with the same parameters that 
were used in the previous calculations. We see that the 
position-momentum uncertainty product decreases and 
flattens out after about 50 measurements. It does not 
reach the minimum value allowed by the Heisenberg un- 
certainty relation but is only about a factor 5-10 times 
larger than this minimum. Examining a little closer the 
individual realizations, we find variations, and plotting 
the 50 percent quantile, i.e., the value which is above 
the uncertainty product of the 50 percent of the realiza- 
tions with the smallest uncertainty products, we observe 
results close to the Heisenberg limit. In these cases we 
are very close to minimum uncertainty states, and we 
know that they are very close to either coherent states 
or squeezed states of motion. With the assumption that 
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FIG. 3: (Color online) The upper curve shows the mean value 
of the position-momentum uncertainty product as a function 
of time obtained from sampling 2000 independent iterations. 
50 % of the realizations have their uncertainty product below 
the lower curve in the figure (the 50 % quantile). The same 
parameters were used for the simulations as in Fig. [5] 

the majority of realizations converge to coherent states, 
we make a simple coherent state fidelity analysis by com- 
puting after 150 measurements for every realization the 
state overlap (a|p|a) between the conditioned state p and 
a coherent state of motion |a), having the same mean 
position and momentum. A histogram of these results 
is shown in Fig. Q] for 500 runs. As can be seen a large 
fraction of the runs results in states which have a co- 
herent state fidelity above 0.9, but we also find a siz- 
able fraction with a nearly vanishing coherent state over- 
lap. To explain these findings, in Fig. [5] the Q-function, 
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FIG. 4: (Color online) Histogram of values of the coherent 
state fidelity after 150 measurements for 500 simulations. 



Q(a) — (a\p\a), of four final states are plotted as func- 
tion of the variable coherent state argument a. Two of 
the examples (a,b) have a high coherent state fidelity and 
naturally show a Q-function that resembles that of a co- 
herent state. The other two (c,d) have a coherent state 
fidelity lower than 0.1. In these cases, we see that the 
main peaks are very similar to the coherent state peaks 
in the other figures. We also see, however, that there are 
areas in the complex plane away from the main peaks 
where the Q-functions are not negligible. The population 
in these regions shift the position and momentum expec- 
tation values so that Q-function may effectively become 
very small at these mean values. The maximum value of 
the Q-function, max a (a\p\a), over coherent state argu- 
ments would thus constitute a more fair measure of the 
coherent state character of the states obtained. 

Our simulations indicate that the measurement back- 
action generates states close to coherent states. To un- 
derstand how this can result from measurements of the 
position only, it is important to recall, that information 
about the position at one moment is converted into infor- 
mation about the momentum at later times due to the dy- 
namics of the system. This is in competition with the ad- 
verse effect that position measurements at the same time 
cause an uncertainty in the momentum. Our measure- 
ments via the Cooper-pair box do not, however, cause 
an instant projection of the oscillator on position eigen- 
states, but rather constitute weak position measurements 
and we thus gradually obtain information about both po- 
sition and momentum of the oscillator. 

To secure a symmetric probing of both position 
and momentum, we have simulated an interaction and 
detection protocol incorporating a wait time of 0.25/oj 
between all measurement. This wait time ensures an 
alternating measurement sequence between the oscil- 
lator quadratures corresponding to equivalent probing 
of position and momentum properties. In Fig. [6] a 
histogram of 500 attempts with this wait time between 
all meauserements is shown. For all 500 attempts the 
coherent state fidelity, defined via the state with the 
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FIG. 5: (Color online) Q-functions for four different simula- 
tions. In (a) and (b) the states are close to coherent states. 
The uncertainty products are 1.2 and 1.0, respectively, in 
units of h/2 and the coherent state fidelities are 0.96 and 
0.99. In (c) and (d) the uncertainty products are 178 and 
109, respectively, in units of h/2, and we have state fidelities 
below 0.1 of coherent states with the same mean position and 
momentum. 



V. FEEDBACK COOLING TO OSCILLATOR 
GROUND STATE 

Since we produce coherent states of the oscillator by 
our measurement back action, and the measurement 
record is known to the experimentalist, it is possible 
to convert the final state into the cantilever motional 
ground-state by application of a force. Since, however, 
our knowledge may not be fully accurate, and the force 
may be imprecise, it may be better to apply a feedback 
after each measurement and to adaptively let the dy- 
namics produce and simultaneously verify the oscillator 
ground state. The feedback force can for example be ac- 
comodated by a piezo-resistive mount for the cantilever 
or by light pressure from a laser beam. It will be repre- 
sented in the Hamiltonian by a term rjx. Since the charge 
state of the Cooper-pair box is known when the feedback 
is applied the interaction with the Cooper-pair box can 
be readily compensated in the feedback term. If the ex- 
pectation values of the position and momentum after a 
measurement procedure are estimated to be (x(0)) and 
(p(0)), the strength -q and the duration r of the feedback 
force are given by, 

(p(0)) 2 +m 2 ^(x(0)) 2 
V 2mco(x(0)) 1 ' 



same mean position and momentum, is now above 0.91, 
as indicated in Fig. [5] 



tan(wr) 



FIG. 6: (Color online) Histogram of the coherent state fideli- 
ties for 500 simulations after 150 measurements have taken 
place. A wait time of 0.25/u; is used between measurements. 



In fact, not only do these wait times give a better re- 
sult, they also represent a more realistic model of the 
measurement procedure allowing a finite read-out time 
for the Cooper-pair box. 



2mu(x(0))(p(0)) 
m 2 uj 2 (x(0)) 2 - (p(0)) s 



(8) 



where m is the effective mass of the oscillator. 

We have simulated this evolution, and in Fig. [7] a his- 
togram of the ground-state population for a hundred it- 
erations is shown. 160 measurements were performed 
in each iteration, and the same system parameters were 
used as for the other simulations. We see that these sim- 
ulations results in a very high ground state population. 



A. The environment 

In the previous sections we treated the oscillator and 
Copper-pair box as if they were perfectly isolated from 
the environment between the measurement processes. 
The cantilever, however, is mechanically mounted, and 
it thermalizes to its environment at a rate determined by 
material and design parameters. The environment will 
thus limit our cooling as it tries to restore thermal equi- 
librium. The coherence time of the Cooper-pair box can 
be made as long as a few fxs [Hj]. Since the Cooper- 
pair box is put into a pure state each time its charge 
state is measured and the time between subsequent mea- 
surements is much smaller than the coherence time, the 
effect of the environments effect on the Cooper-pair box 
can be neglected. Therefore we concentrate on environ- 
mental effects on the cantilever mode, which we model 
by damping terms in a master equation for the oscillator 
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FIG. 7: (Color online) The ground-state fidelity after 160 
measurement and feedback operations. The simulations were 
performed with the same parameters as the ones shown in 
Fig. ® 



FIG. 8: (Color online) Steady state excitation of the cantilever 
mode as a function of the thermal occupation of the reservoir 
for three different Q-values. 



density matrix. If H is the system Hamiltonian, the ap- 
propriate master equation describing coupling with rate 
7 to a black-body reservoir with a mean excitation of n 
quanta at the cantilever resonance frequency, is given by 

dp — % 7 

— = — [H,p] + -(n+l)(2apa [ - a) ap - pa) a) 

+—n(2a^ pa — aa) p — pda)) (9) 

From this equation we expect that both the thermal oc- 
cupation of the reservoir n and the coupling constant 7 
will play a role in the dynamics of the oscillator. We 
make the assumption that after a 100 measurements the 
system has reached steady state. Note that this is not a 
stationary state of the system, since each measurement 
and feedback process may both cause net excitation and 
de-excitation of the oscillator. In steady state, the popu- 
lation of the oscillator eigenstates after each feedback will 
thus fluctuate around a certain " typical" mean value. We 
expect that the memory of the oscillator of these steady 
state fluctuations will be limited to a couple of rounds 
of measurements and feedback, and hence the mean en- 
ergy in steady state is effectively sampled by a long time 
average. 

Instead of the coupling strength it is customary to 
specify the Q-value for oscillators. This is related to the 
coupling constant by Q — ui/2-f. Micromechanical oscil- 
lators have typical Q-values between 10 3 and 10 5 . In Fig. 
[S] we show the mean steady state energy of the oscillator 
plotted against the thermal occupation of the reservoir 
for three different Q-values. Fig. [8] provides a measure 
of how well the cooling counters the environmental heat- 
ing. For a Q-factor of 10 5 we get something very close 
to ground-state cooling. The mean excitation (n) is al- 
most zero so that the occupation of the ground state is 
almost unity and the results show that the environment 



a mean excitation well below unity, while for Q = 10 3 
the mean occupation is on the order of unity but still far 
below n of the reservoir. 



VI. CONCLUSION 

In summary, we have proposed to extract energy from 
a micromechanical oscillator by read-out of its position 
by means of a Cooper-pair box, and by acting back on 
the state of the oscillator according to the read-out. Un- 
like cooling schemes which explicitly damp the cantilever 
by the coupling to the readout device, our basic interac- 
tion does not have any such average cooling effect. The 
mechanism in fact extracts entropy rather than energy, in 
the sense that the interaction and measurements extract 
information about the cantilever and hence reduces its 
entropy until we reach a pure quantum state. By suit- 
able timing of the measurements, this pure state hap- 
pens to be very close to a coherent state, whose mean 
position and momentum are revealed by the detection 
record, and which can hence be deterministically dis- 
placed to the ground state - either at the end of the 
measurement record or by small steps after each mea- 
surement event. In comparison to the cooling schemes 
that work by a genuine damping mechanism, our cooling 
may work for a wider range of frequency and coupling 
parameters, and even for parameters where the average 
back action effect causes heating. Provided, of course, 
that we can provide the feedback force on the cantilever. 
In [l[ a scheme is presented where similar physical mea- 
surements can be used to generate non-classical states of 
an oscillator which is initially in a coherent state. Our 
analysis uses the same architecture, and hence we pro- 
vide a means to obtain the initial state needed for that 
and other protocols for generation of non-classical states. 
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